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Abstract. Results of experiments investigating air showers in the energy region of the knee 
are summarized. The all-particle energy spectrum, the mean logarithmic mass, and the average 
depth of the shower maximum will be discussed. Spectra for groups of elements from air shower 
' data are compared to results from direct measurements. 
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' 1. Introduction 

QQ . One of the most remarkable structures in the energy spectrum of cosmic rays is a change of 

I the spectral index 7 of the power law dN/dE oc E'^' at an energy of about 4 PeV, the so called 

in . knee. In the literature various reasons for the origin of the knee are discussed Some popular 

\ concepts are: The maximum energy attained during the acceleration at supernova remnant shock 

' waves 01 m EI El ) reacceleration in the galactic wind [7] , acceleration in pulsars |H1 , or leakage 

I . of particles from the Galaxy during the propagation process (HI El El El E] • Different 

O \ scenarios are the acceleration of cosmic rays in 7-ray bursts jl5l 1161 117j . During the propagation 

' through the Galaxy cosmic-ray particles are proposed to interact with dense photon fields close 

d . to the sources or with background neutrinos 1181 1191 12UI 121 1 . As a last possibility, new processes 

between elementary particles in the atmosphere are postulated which may transfer energy into 
non-observed channels and thus cause a knee in the observed air shower components |221 1231 124j . 
An answer to the question of the cause for the icnee is expected to reveal also crucial 
' information on the origin of galactic cosmic rays. Experimental access to such questions is 

provided by measurements of charged cosmic rays (the classical nucleonic component) and 7- 
rays with experiments above the atmosphere and by the observation of air showers initiated by 
high-energy particles in the atmosphere. 

A wealth of information on potential cosmic-ray sources is provided by recent measurements 
of TeV 7-rays from shell type supernova remnants with the H.E.S.S. experiment j25j . The 
observations reveal a shell structure of the remnant and an energy spectrum of 7-rays oc E""^'^ 
in agreement with the idea of particle acceleration in the shock front. The spectrum extends up 
to energies of 10 TeV and provides evidence for the existence of particles with energies beyond 
100 TeV at the shock front that emerged from the supernova explosion. 

The gyromagnetic radius of a proton with an energy of 1 PeV in the galactic magnetic fields is 
about 0.4 pc and consequently the probability to detect a source in charged galactic cosmic rays 
is rather small. Indeed, observations by the KASCADE experiment did not reveal any significant 

^ Invited talk, presented at the "Workshop on Physics of the End of the Galactic Cosmic Ray Spectrum", Aspen, 
USA, April 25 - 29, 2005 
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Figure 1. All-particle energy spectrum. Results from direct measurements by Grigorov et al 
inni, JACEE |HI], RUNJOB |S2], and SOKOL [HHl as well as from the air shower experiments 
AGASA [H3 Akeno 1 km^ ^ and 20 km^ [HE], AUGER BASJE-MAS [HHl, BLANCA 
|Sn], CASA-MIA 001, DICE HI], EAS-TOP Fly's Eye US], Haverah Park (1991) ^ and 
(2003) HE], HEGRA jMl, HiRes-MIA [H], HiRes-I HE], HiRes-II HH], KASCADE eiectrons 
and muons interpreted with two hadronic interaction models jEO]) hadrons [EJ, and a neural 
network analysis combining different shower components jHS], MSU jSH]; Mt. Norikura |54j . 
SUGAR Tibet AS7 [HE] and AS7-III [HZ], Tunka-25 p], and Yakutsk The hues 

indicate the spectrum according to the poly-gonato model. 



evidence for point sources in the energy range from 0.3 to about 100 PeV (2^]. Nevertheless, 
information on the sources is derived from detailed measurements of the cosmic-ray composition. 
Investigations of the abundance of refractory nuclides reveal that their abundance at the sources 
is extremely similar to the abundance observed in the solar system ji27j . This indicates that 
cosmic rays are accelerated from a well- mixed sample of contemporary interstellar matter. 

Information on the propagation of particles through the Galaxy is obtained by measurements 
of the ratio from primary to secondary nuclei, the latter are produced in spallation reactions of 
cosmic rays with particles of the interstellar medium, and the abundance of radioactive nuclides 
in cosmic rays. Measurements by the CRIS experiment yield a residence time for cosmic rays in 
the Galaxy of 15-10^ yrs and a propagation path length of about 10 g/cuP for particles with GeV 
energies |28^. The measured abundances at GeV energies are frequently described using leaky- 
box models. At higher energies measurements of anisotropy amplitudes yield information on 
the propagation process [IHl- They indicate that the propagation can be described by diffusion 
models, while leaky-box models exhibit too large anisotropy values even at relatively low energies. 

At PeV and EeV energies, information is derived from the measurements of the energy 
spectrum and mass composition of cosmic rays. In the following, experimental results are 
reviewed on the all-particle energy spectrum (Sect. {2]), the average mass composition (Sect.jHJ), 
and spectra for elemental groups (Sect.^. 

2. All-particle energy spectrum 

The all-particle energy spectra obtained by many experiments are compiled in Fig. ^ Shown 
are results from direct measurements above the atmosphere as well as from various air shower 
experiments. The individual measurements agree within a factor of two in the flux values 
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Figure 2. Mean logarithmic mass of cosmic rays derived from the measurements of electrons, 
muons, and hadrons at ground level. Results are shown from CASA-MIA jHJ, Chacaltaya |62j . 
EAS-TOP electrons and GeV muons [fil, EAS-TOP/MACRO (TeV muons) [Ml, HEGRA CRT 
|65j ■ KASCADE electrons and muons interpreted with two hadronic interaction models [BOj . 
hadrons and muons [MI; as well as an analysis combining different observables with a neural 
network .52 , and SPASE/AMANDA [HTj. The lines indicate expectations according to the 
poly-gonato model. For comparison, results from direct measurements are shown as well from 
the JACEE [HH] and RUNJOB ^ experiments. 



and a similar shape can be recognized for all experiments with a knee at energies of about 
4 PeV. Typical values for the systematic uncertainties of the absolute energy scale for air shower 
experiments are about 15 to 20%. Renormalizing the energy scales of the individual experiments 
to match the all-particle spectrum obtained by direct measurements in the energy region up 
to almost a PeV requires correction factors in the order of ±10% [Sni- A remarkable result, 
indicating that behind an absorber of 11 hadronic interaction lengths or 30 radiation lengths 
the energy of the primary particle is determined with an absolute error in the order of ±10%. 
One should keep in mind that the experiments investigate different air shower components, are 
situated at different atmospheric depths, and use different interactions models to interpret the 
observed data. Nevertheless, the systematic differences are relatively small and the all-particle 
spectrum seems to be well known. Average experimental fluxes are tabulated in |60j . 

The lines in Fig. Q indicate sum spectra obtained by extrapolations of the energy spectra for 
individual elements from direct measurements assuming power laws with a cut-off proportional 
to the charge of the respective element according to the poly-gonato model ||6Dj. Sum spectra for 
elements from hydrogen to nickel and for all elements are shown. The extrapolated all-particle 
flux is compatible with the flux derived from air shower experiments in the knee region. Above 
10^ GeV the flux of galactic cosmic rays is not sufficient to account for the observed all particle 
spectrum, and an additional, presumably extragalactic component is required. 

3. Average mass of cosmic rays 

At energies below a PeV energy spectra for individual elements have been observed above the 
atmosphere [Hni- At higher energies this is presently not possible due to the low flux values and 
the large fluctuations in the development of extensive air showers. Thus, mostly the mean mass 
is investigated. An often-used quantity to characterize the composition above 1 PeV is the mean 
logarithmic mass, defined as {In A) = J2i i^i lii^j, being the relative fraction of nuclei of mass 
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Figure 3. Average depth of the shower maximum X^ax as function of primary energy as 
obtained by BLANCA (HHI, CACTI [Zl], DICE gl], Fly's Eye Haverah Park [72 , HEGRA 
HE], HiRes/MIA j2B], Mt. Lian Wang SPASE/VULCAN HHj, Tunka-25 and Yakutsk 
|76j . The Unes indicate simulations for proton and iron induced showers using the CORSIKA 
code with the hadronic interaction model QGSJET 01 (solid line) and a version with lower cross 
sections and slightly increased elasticities (dashed line, model 3 in |77j ) . 

Ai. In the superposition model of air showers, the shower development of nuclei with mass A 
and energy Eq is described by the sum of A proton showers of energy E = Eq/A. The shower 
maximum t penetrates into the atmosphere as t oc Inii^, hence, most air shower observables at 
ground level scale proportional to In A. 

Frequently, the ratio of the number of electrons and muons is used to determine the mass 
composition. The experiments CASA-MIA (HH, EAS-TOP or KASCADE use muons with 
an energy of several 100 MeV to 1 GeV. To study systematic effects two hadronic interaction 
models are used to interpret the data measured with KASCADE [SIl- High energy muons 
detected deep below rock or antarctic ice are utilized by the EAS-TOP/MACRO [H^ and 
SPASE/ AMANDA [HZI experiments. Also the correlation between the hadronic and muonic 
shower components has been investigated by KASCADE 66 . The production height of muons 
has been reconstructed by HEGRA/CRT |M] and KASCADE jZO]. 

Results from various experiments measuring electrons, muons, and hadrons at ground level 
are compiled in Fig. [2l At low energies the values for the mean logarithmic mass are compared 
to results from direct measurements. A clear increase as function of energy can be recognized. 
The experimental values follow a trend predicted by the poly-gonato model as indicated by the 
lines in the figure. However, individual experiments exhibit systematic deviations of about ±1 
unit in {In A) from the line. Of particular interest are also the investigations of the KASCADE 
experiment, interpreting the same measured data with two different models for the interactions 
in the atmosphere results in a systematic difference of about 0.7 to 1 in (lnj4). 

Another technique to determine the mass of cosmic rays are measurements of the average 
depth of the shower maximum using non-imaging and imaging Cerenkov detectors as well as 
fluorescence telescopes at the highest energies. The results of several experiments are presented 
in Fig. 121 as function of energy. The observed values are compared to predictions of air shower 
simulations for primary protons and iron nuclei using the program CORSIKA ^78_ with the 
hadronic interaction model QGSJET 01 ^79^ and a modified version with lower cross sections 
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Figure 4. Mean logarithmic mass of cosmic rays derived from the average depth of the shower 
maximum, see Fig. |21 Two hadronic interaction models are used to interpret the measurements: 
QGSJET 01 (top) and a modified version with lower cross sections and a slightly increased 
elasticity (model 3a 77 , bottom). For references, see caption of Fig. |3J For comparison, 
results from direct measurements are shown as well from the JACEE [HSI and RUNJOB |32j 
experiments. The lines indicate expectations according to the poly-gonato model. 



and larger values for the elasticity of the hadronic interactions (model 3a in [77]). The latter is 
compatible with measurements at colliders, for details see jTTj. The lower values for the total 
inelastic proton-antiproton cross sections yield also lower inelastic proton-air cross sections, 
which are in good agreement with recent measurements from the HiRes experiment I81j. 
In principle, the difference between the solid and the dashed lines in the figure represents an 
estimate of the projection of the experimental errors from collider experiments on the average 
depth of the shower maximum in air showers. At 10^ GeV the difference between the two 
model versions for primary protons is about half the difference between proton and iron induced 
showers. This illustrates the significance of the uncertainties of the collider measurements for 
air shower observables. 

Knowing the average depth of the shower maximum for protons X^^^j. and iron nuclei X^'^^^. 
from simulations, the mean logarithmic mass is derived in the superposition model of air showers 
from the measured XZaT using (In^) = {XZaT " XP,,,)/iX^l, - X^^J • InAp^. The 
corresponding {In A) values, obtained from the results shown in Fig. |21 are plotted in Fig. 
versus the primary energy using both interaction models to interpret the observed data. Up to 
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Figure 5. Energy spectrum for protons. Results from direct measurements above the 
atmosphere by AMS ^22,, ATIC 84 , BESS CAPRICE |HE], HEAT [HZI, Ichimura et al 
(HE], IMAX JACEE [nni, MASS [HI], Papini et al RUNJOB |S2j, RICH-II [HHl, Ryan 
et al |94j . Smith et al '95', SOKOL "33^, Zatsepin et al |96j . and fluxes obtained from indirect 
measurements by KASCADE electrons and muons for two hadronic interaction models |5f)j and 
single hadrons [HZI, EAS-TOP (electrons and muons) [HHI and single hadrons [22], HEGRA |innj . 
Mt. Chacaltaya ToT, Mts. Fuji and Kanbala HOSj, Tibet burst detector (HD) and AS7 

(HD) |ilfl4|| . The line indicates the spectrum according to the poly-gonato model. 

about a PeV there are only marginal differences between the two interpretations. On the other 
hand, at large energies a significantly heavier composition is obtained for the modified version 
(model 3a). At 10^ GeV the differences amount to about 1 in (In A). These examples illustrate 
how strong the interpretation of air shower measurements depends on model parameters such 
as the inelastic cross sections or elasticities used. At Tevatron energies the cross sections vary 
within the error range given by the experiments and at 10*^ GeV the proton air cross sections 
of QGSJET and model 3a differ only by about 10%, but the general trend of the emerging 
{In A) distributions proves to be significantly different. This underlines the importance to test 
and improve the understanding of hadronic interactions in the atmosphere with air shower 
experiments |H2l- R seems that the interpretation with model 3a is better compatible with the 
mean logarithmic masses derived from electrons, muons, and hadrons already shown in Fig. [2l 
and also with the predictions of the poly-gonato model as indicated by the lines. 

4. Spectra for elemental groups 

A significant step forward in the understanding of the origin of cosmic rays are measurements 
of energy spectra for individual elements or at least groups of elements. Up to about a PeV 
direct measurements have been performed with instruments above the atmosphere. As examples, 
results for primary protons, helium, and iron nuclei are compiled in Figs. |Slto[7| Recently, also 
indirect measurements of elemental groups became possible. 

A special class of events, the unaccompanied hadrons were investigated by the EAS-TOP and 
KASCADE experiments |991l97j . Simulations reveal that these events, where only one hadron is 
registered in a large calorimeter, are sensitive to the flux of primary protons. The derived proton 
fluxes agree with the results of direct measurements as can be inferred from Fig. El indicating 
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Figure 6. Energy spectrum for helium nuclei. Results from direct measurements above the 
atmosphere by ATIC [H3, BESS 85 , CAPRICE |HE], HEAT |H3, Ichimura et al [HHI, IMAX 
jHi, JACEE [201, MASS inU, Papini et al RICH-II (M!, RUNJOB |S2], Smith et al [Ml, 
SOKOL I2HI) Webber |in5j . and fluxes obtained from indirect measurements by KASCADE 
electrons and muons for two hadronic interaction models Mts. Fuji and Kanbala 102 , and 
Tibet burst detector (HD) |in3j . The line indicates the spectrum according to the poly-gonato 
model. 



a reasonably good understanding of the hadronic interactions in the atmosphere for energies 
below 1 PeV. 

At higher energies a breakthrough has been achieved by the KASCADE experiment. 
Measuring simultaneously the electromagnetic and muonic component of air showers and 
unfolding the two dimensional shower size distributions, the energy spectra of five elemental 
groups have been derived [SOI- In order to estimate the influence of the hadronic interaction 
models used in the simulations, two models, namely QGSJET 01 and SIBYLL |in6j . have been 
applied to interpret the measured data. It turns out that the all-particle spectra obtained agree 
satisfactory well within the statistical errors. For both interpretations the flux of light elements 
exhibits individual knees. The absolute flux values differ by about a factor of two or three 
between the different interpretations. However, it is evident that the knee in the all-particle 
spectrum is caused by a depression of the flux of light elements. The KASCADE results are 
illustrated in Figs. El to [Tj 

In the figures also results from other air shower experiments are shown. EAS-TOP derived 
spectra from the simultaneous observation of the electromagnetic and muonic components. 
HEGRA used an imaging Cerenkov telescope system to derive the primary proton flux [lOOj . 
Spectra for protons and helium nuclei are obtained from emulsion chambers exposed at Mts. 
Fuji and Kanbala jl02j . The Tibet group performs measurements with a burst detector as well 
as with emulsion chambers and an air shower array jlU3l I57j . 

Over the wide energy range depicted, the measurements seem to follow power laws with a 
cut-off at high energies. The spectra according to the poly-gonato model are indicated in the 
figures as lines. It can be recognized that the measured values are compatible with cut-offs at 
energies proportional to the nuclear charge Ez = Z ■ 4.5 PeV. 
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Figure 7. Energy spectrum for iron nuclei. Results from direct measurements above the 
atmosphere by CRN [M], HEAO-3 [MI, Juliusson et al [THi, Minagawa et al [1101, TRACER 
(single element resolution) and Hareyama et al |112j . Ichimura et al (HE]) JACEE pT] . 
RUNJOB [21], SOKOL (iron group), as well as fluxes from indirect measurements (iron 
group) by KASCADE electrons and muons for two hadronic interaction models |5Uj and EAS- 
TOP The line indicates the spectrum according to the poly-gonato model. 



5. Conclusion and outlook 

In the last decade substantial progress has been achieved and the knowledge about galactic 
cosmic rays has been significantly increased. The all-particle energy spectrum is reasonably 
well known. For the first time, spectra for groups of elements could be derived from air shower 
measurements. The observed spectra seem to exhibit cut-offs proportional to the nuclear charge. 
The increase of the mean logarithmic mass derived from air shower observations seems to be 
compatible with subsequential cut-offs for individual elements. However, the interpretation of 
the measurements is still limited by the uncertainties of the description of hadronic interactions 
in the atmosphere. The implications of recent measurements on the contemporary understanding 
of the origin of the knee have been discussed elsewhere |113j . 

With the recent measurements of TeV 7-rays from supernova remnants a new window 
has been opened, and the hypothesis of cosmic-ray acceleration in supernova remnants is tested 
directly. 

With the KASCADE-Grande experiment |114j . investigations of the energy spectra for groups 
of elements will be extended into the region of the expected iron icnee (~ 100 PeV) and up to 
the second knee (~ 400 PeV). In particular, the inclusion of the hadronic component will help 
to improve the interaction models at these energies. 

A new technique to measure air showers is about to be established: The LOPES experiment 
has registered first air showers observing geosynchrotron emission in the radio frequency range 
from 40 to 80 MHz |TT^ . 

The presently largest balloon borne detector, the TRACER experiment almost reaches 

the energy region of indirect measurements. Successor experiments will provide twofold progress 
for the understanding of the origin of galactic cosmic rays: an extension of the measurements of 
the ratio of secondary to primary nuclei to energies approaching the knee will improve the 
knowledge about the propagation of cosmic rays and the extension of energy spectra with 
individual element resolution towards the air shower regime will provide useful information 



to improve hadronic interaction models. 
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